Density of state and non-magnetic impurity effects in electron-doped cuprates 
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we analyze the density of state (DOS) and a non-magnetic impurity effect in electron-doped 
cuprates starting from two different scenarios: the d^2_y2-wave superconductivity coexisting with 
antiferromagnetic spin density wave (SDW) order versus d,^2_y2-wa,ve superconductivity with a 
higher harmonic. We find that in both cases the local density of state (LDOS) exhibits two impurity- 
induced resonance states at low energies. We also find that for the intermediate value of the SDW 
gap, the DOS looks similar to that obtained from the scenario of the d^2_y2-wave gap with a higher 
harmonic, suggesting the presence of a non-monotonic d^2_j^2-wave gap. However, if the SDW gap 
is sufficiently large the DOS looks more conventional s-wave like. This obvious difference from the 
DOS resulted from the d^2_y2-wave gap with a higher harmonic model, could differentiate the two 
above scenarios and is needed to be proved in the further doping dependence of tunneling spectrum 
measurement. 

PACS numbers: 74.72.Jt, 74.20.Mii, 74.25.Ha, 74.25.Jb 
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I. INTRODUCTION 

The electron-doped cuprate high-temperature super- 
conductors have recently attracted a considerable at- 
tention due to their considerable asymmetry from hole- 
doped ones. For instance, it has been known that in con- 
trast to the hole-doped cuprates, the electron-doped sys- 
tems show (i) relatively low superconducting transition 
temperature, (ii) narrow superconducting (SC) region in 
the phase diagram, and (iii) the robust antiferromagnetic 
(AF) orderi 

The phase sensitive scanning SQUID measurements^, 
nuclear magnetic resonance study-, and also ARPES 
experiments^t^ have provided strong evidences that the 
electron-doped superconductors are the dx2-y2-wa.ve su- 
perconductors. This conclusion seems to be natural, 
since from the theoretical point of view, the same pairing 
mechanism is expected both for hole- and electron-doped 
materials, although the functional form of the dx2-y2- 
wave gap in electron-doped materials is more subtle is- 
sue. The high-resolution ARPES data on the leading- 
edge gap in PrQ,sgLaCeo,iiCu04'^ , Raman experiments 
in NCCO-, and doping dependence of tunneling study 
in Pr2-xCexCuOi-s^ show a non-monotonic dj.2_y2- 
wave gap with a maximum value in between nodal and 
antinodal points on the Fermi surface (FS), and mean- 
while the measurements of optical conductivity cri[uj) in 
P'ri.s5CGo.i5CuO/^ were also interpreted as an indirect 
evidence of a non-monotonic dx2-y2-wsLve gap. 

The theoretical explanations for the origin of the non- 
monotonic gap behavior are classified into two scenar- 
ios so far. One is the coexisting AF with SC order 
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the non-monotonic d^2_j^2-wave gap appears naturally 
under the assumption that the dx2_y2--wsLve pairing is 
caused by the interaction with the continuum of over- 
damped AF spin fluctuations. Spin-mediated interaction 
is attractive in the dx2-y2-wa,ye channel and yields a gap 
which is maximal near the hot spots - the points along 
the FS, separated by AF moment Qaf- In optimally 
doped NCCO and PCCO, hot spots are located close to 
Brillouin zone diagonals, and one should generally expect 
;-wave gap to be non-monotonic. In this case. 



the d^: 

the non-monotonic gap behavior is an intrinsic property 
in the SC state regardless of the presence of the AF order. 

In fact, impurity effect has always been used as one 
of the most important and effective tools to distinguish 
pairing symmetry in the conventional and unconventional 
superconductors. It has been known that the Yu-Shiba- 
Rusinov state in the conventional BCS superconductor 
(s-wave) was located at the gap edge^'^. However in the 
d,j.2_y2--wave superconductor an impurity-induced bound 
state appear near the Fermi energyl^. The origin of this 
difference has been explained as different phase structure 
of two kinds of SC order: d^2_y2-'wave pairing symme- 
try with the line nodal gap while s-wave symmetry with 
nodeless gap. Therefore, we in this paper analyze the 
DOS and a non-magnetic impurity effect starting from 
two different scenarios: the (i2,2_j^2-wave superconductiv- 
ity coexisting with AF SDW wave order versus dx2_y2- 
wave superconductivity with a higher harmonic. We 
find that in both cases the LDOS exhibits two impurity- 
induced resonance states at low energies. We also find 
that for the intermediate value of the SDW gap, the 
DOS looks similar to that obtained from the scenario 
scenaric^'^". Although the SC gap itself is assumed to of the d^^ -y2-wa.ve gap with a higher harmonic, suggest- 



be monotonic d^2_y2-wave, when the AF order is in- 
troduced, the resulting quasiparticle excitation can be 
gapped by both orders and behaves to be non-monotonic 
(i2,2_y2-wave gap. The other one is to extend the SC 
gap out of the simplest dx2_y2-waYe via the inclusion of 
a higher harmonic termiiii^. From theory perspective. 



ing the presence of a non-monotonic d^2_y2-wave gap. 
However, if the SDW gap is sufficiently large the DOS 
looks more conventional s-wave like instead of the non- 
monotonic dx2-y2-wave gap behavior. This obvious dif- 
ference from the DOS resulted from the dx2_y2-wa.ve gap 
with a higher harmonic model even in the underdoped 
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regimes, could differentiate the two above scenarios. Our 
results strongly suggest the further doping dependence 
of tunneling spectrum measurement, especially in the 
undcrdoped regimes, is needed to be carried out so as 
to shed light on the physical origin of the unusual non- 
monotonic gap in the electron-doped cuprates. 



II. THE MODEL AND T-MATRIX 
FORMULATION 

We start from a phenomenological superconducting 
Hamiltonian on a square lattice at the mean field level, 



where 



H 

Hsc 



Hsc + HsDWi 
^[(ek-M)cLck<, 



kcr 



Hs 



DW 



^k(4TC-ki +C-kiCkT)], 
- ^ Mcr(c[^Ck+QCT + /i.e.). 



(1) 



where c\^^ (cko-) is the fermion creation (destruction) op- 
erator and /i is chemical potential. The SC order param- 
eter has the usual form Ak = A(cos(fca;) — cos(fcj,))/2 and 
the AF SDW order parameter is M . We take the normal 
state tight binding dispersion ek, which is written as 

Ek = — 2t(cOs(fcj,) -f- COs(fcj,)) — 4^1 COs(fca:) COs(fcj^) 

— 2t2(cos(2fca;) -f cos(2fcj,)) 

— 4f3(cos(2fca;) cos(fcy) -I- cos(fc2,) cos(2fcj^)) 

— 4t4COs(2fca;) cos(2fc2;) (2) 

with t = 120 meV, h = -60 meV, ^2 = 34 meV, h = 7 
meV, ti — 20 meV, and /i = -82 meV at 0.11 dopingl^. 
reproducing the underlying FS as inferred from recent 
ARPES experiment^. Note that here the wave vector k 
is restricted to the magnetic Brillouin zone (MBZ). 

It is convenient to introduce a 4 x 4 matrix for- 
mulation to treat the coexisting SDW and SC phase. 
After taking a four-component Nanibu spinor ^'k = 

)^ with Q = (tt, tt) being the 



-k-Qi 



(ckT,Ck+QT,c^ 

nesting vector, we rewrite the Hamiltonian as 
H = Y^ *^ ((ek - m)t-3/Oo + MtiPo + Akr3pi)*k, (3) 



where T^pi 





T3 



T3 





. Then the single-particle matrix 



Green's function is determined as Gq (fc, iujn) = iuJn — 
(ffk - m)t3Po - MtiPo - AkTspi. 

The scattering of quasiparticles from the impurity is 
described by a T-matrij*i^ii^iii, T{iujn), which is local 
and independent of wave vectors. Thus, we define the 
2x2 Green's function in the presence of a single impurity 
as 



G{i,j;iujn) = Co(i -i;iw„) 



(4) 



Co(«, j;«^^«) 



with 



N ^ 

k 

Gj(fc,iw„) 
Gg(fc,iw„) 



e'Q-^'J[Go]f (A:,zc.„; 



G^(fc,itj„) 



(5) 



[Golf (A:,zt^„) 



[Go]i^(fc 



(6) 



Here I — 1,2,3,4 denotes the left-top, right-top, 
left-bottom and right-bottom 2x2 block element of 
Go{k,iujn), and Ri is lattice vector and Ry = Ri — Rj. 
The T-matrix can be written by 



T{iuj„) 



UoP3 



1 - UoP3Co{0,0;iu>n) 



(7) 



For the d-wave pairing symmetry, one can find that the 
local Green's function Co(i,«;iw„) is diagonal. As a re- 
sult, the diagonal T-matrix 



ril,22(»t^n) 



±Uo 



l-C/o[Co(0,0;±za;„)] 
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(8) 



where the upper (lower) sign denotes Tn (122), will give 
rise to a particle- (cures < 0) and hole-like (cj^es > 0) 
resonance state. These resonance states generate the 
sharp peaks in the LDOS only in the unitary limit 

(I LUres | /A < 1) whcrC 1 = UoRe[Co{0 , 0; ±UJres)]ll ■ 

Finally, these above equations allow a complete solution 
of the problem as long as the order-parameter relaxation 
can be neglected. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

A. d^2_y2-wave coexisting with AF SDW order 



We firstly calculate the DOS 

p{uj) = --Imy^ Gii{k,uj) 



(9) 



plotted in Fig.l for different SDW gap (M). In princi- 
ple, SDW gap (M) and SC gap (A) need to be solved 
self-consistently. For simplicity, we in the following cal- 
culation choose A — O.OOSey reproducing the ARPES 
estimate of SC gap in Pro,sQLaCeo,iiCu04^, and dif- 
ferent values of SDW gap M = 0.14,0.12,0.05, and 
0{eV), where the maximal value M — 0.14(ey) is a self- 
consistent value from Ref.15 at doping x — 0.11, and 
the gradually decreasing values with the independent 
band dispersion (Eq.(2)) corresponds to the increasing 
dopingi^iiSiii^. In this case, our SDW gap dependence of 
DOS can be qualitatively equivalent to the doping evo- 
lution of DOS. 
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In the coexisting AF SDW and SC state, there are 
three features of interest. First, at low energies the DOS 
spectrum remains d-wave hke superconductor behavior 
for fairly large values M . Second, with the increasing 
of M, besides the coherent peak located at the SC gap, 
another two new "coherent peaks" appear symmetrically 
at positive and negative low energy, and the positions 
of these peaks shift towards the Fermi energy, indicat- 
ing the presence of non-monotonic da;2_y2-wave gap be- 
havior. The third and last feature is that at sufficiently 
large SDW gap M, the DOS shows the U-shaped behav- 
ior instead of the non-monotonic (ia;2_y2-wave gap be- 
havior. This U-shaped DOS has been observed in earlier 
point contact tunneling spectrum and was explained as 
the character of s-wave pairing^-^"-^-'^, but in present case 
it is the result of the coexisting SDW and SC orders, the 
difference between them will be differentiated by consid- 
ering a non-magnetic impurity effect. 
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FIG. 1: (color online) The DOS for various SDW gap M. 
From bottom to top the SDW gap are successively 0.14, 0.12, 
0.05, and 0(eV). 

These unusual SDW gap dependence of DOS can be 
understood as follows. After diagonalizing the Hamilto- 
nian (1), we obtain the poles of G{k,uj), quasiparticlc 



energy bands i-E^ : 



with 



4" 



Ek + Ek-I-Q _|_ , / (£k - ^k+q) 



Af2 



(10) 



(11) 



In the normal state, the Fermi surface for the corre- 
sponding AF bands (f^ — n) has been plotted in Fig. 2, 
where the dotted-dashed curve is given by the band " -I- " 
and the dotted curve by band " — " . At the intermediate 
value of the SDW gap M = 0.05, the FS has been sepa- 
rated into three parts locating around (tt, 0), (0,7r) and 
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FIG. 2: (color online) Fermi surface for AF bands (^j^ — /i) 
with the SDW gap M = 0.05. The dotted-dashed curves are 
contributed by the band " -I- " and the dotted curve by band 
" — ". The zero-energy contour of £k (solid curves) and £k+q 
(dashed curves) are shown for reference. 



(7r/2, 7r/2) respectively, which is quite agreement with the 
ARPES experiment on NCCO at the optimal doping^^ 
and the previous theoretical calculations^-^-- . In this 
case, the DOS in the SC state obviously shows two dif- 
ferent coherent peaks: one at the low energy corresponds 
to the SC gap maximum opened on the FS piece around 
(7r/2,7r/2), and the other at high energy corresponds to 
FS pieces around (7r,0), and (0, tt). Therefore the non- 
monotonic (ij.2_j,2-wave behavior gap mentioned above 
can be satisfactorily explained to the effective gap in- 
duced by coexisting AF SDW and SC order. With the 
decreasing of the SDW, in the limit case M=0 (solid curve 
in Fig. 2), the two bands in the FS merge into a single one 
again, therefore, the DOS shows a V-shaped behavior, 
similar to that of the hole-doped ones with a monotonic 
gap^i^i^. In contrast, when the SDW increases to be 
sufficiently large (not shown here), the FS piece around 
(7r/2,7r/2) disappears. As a result, the DOS denoted by 
dotted-dashed line in Fig.l behaves to be U-shaped like. 
As expected above, the SDW gap dependence of FS is 
qualitatively similar to the intriguing doping evolution 
of FS^. 

In the presence of a single non-magnetic impurity, the 
LDOS which can be measured in the STM experiment 



TT 

+ -ImG22(r,r;-cj-i0^) (12) 

TT 

with the subscripts 11 and 22 refer to the electron and 
hole part of the diagonal Nambu Green's function, has 
been plotted in Fig.3(a)-3(d) for different AF SDW gap, 
where black solid curve denotes the clean case which is 
equivalent to DOS, and red dashed curve denotes the 
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FIG. 3: (color online)The LDOS for various SDW gap M in 
the clean case and the nearest neighbor of a non-magnetic 
impurity site. (a)M = eV, (b)M = 0.05 eV, (c)M = 0.12 
eV, (d)M = 0.14 eV. 



nearest neighbor of the impurity site. In the calcula- 
tion, we take the same parameters as in Fig.l and set 
the scattering strength [/q — 5eV. Note that our results 
remain qualitatively unchanged when Uq is varied. Since 
the DOS remains c?a;2_y2-wave like at low energies for 
fairly large values M, the two impurity- induced reso- 
nance peaks in LDOS appear at positive and negative 
energy and have different height due to the particle-hole 
asymmetry, similar to the hole-doped cases^^. As M in- 
creases, due to the presence of new "coherence peaks" 
and an inward shift in the position of them in DOS, 
two impurity-induced resonance peaks correspondingly 
shift towards the lower energy and eventually merge into 
one peak at Fermi energy. For the sufficiently large 
M, although the DOS behaves to be U-shaped like, two 
impurity-induced resonance peaks still exist, which is dif- 
ferent from the pure s-wave superconductors where no 
resonance state within the gap was induced by a non- 
magnetic impurity^^ . Thus we come to a conclusion that 
U-shaped DOS is due to the coexisting SDW and SC 
order. We notice that the similar result has also been 
obtained in the recent paper by Wan et al starting with 
t - J model^. 



B. dj.2_j^2-wave with a higher harmonic 

In this scenario, we only take the Hamiltonian Hsc 
without AF SDW order term. Based on a scenario in 
which superconductivity arises from the exchange of AF 
spin fluctuations, it was argued that the maximum SC 
gap is achieved near those momenta on the FS (hot spots) 
that are connected by Q = (tt, 7r)^iiii2&. In the electron- 
doped cuprates, the hot spots are located much closer to 
the zone diagonal, leading to a non-monotonic behavior 



FIG. 4: (color online) (a)The doping (see inset of (b)) depen- 
dence of the FS in the one quarter of Brillouin zone. MBZ is 
denoted by the dashed line and (j) is the FS angle, (b) cor- 
responding non-monotonic gap as a function of the FS angle 
0. 



of the SC gap. A good fit of Ak to the experimental 
data is achieved via the inclusion of a higher harmonic, 
such that Ak — Ao(cos(fca:) — cos(fcy))/2 — Ai(cos(3fcj:) — 
cos(3fcy))/2 with Ao = 5.44 meV and Ai = 2.34 meV en- 
sures that the maximum of |Ak| along the FS is located 
at the hot spots, as shown by the red dashed curve in 
Fig. 4(a) and Fig. 4(b). By tuning the chemical potential 
and Aq, the doping dependence of FS and correspond- 
ingly non-monotonic gap as a function of the FS angle 
(f) have also been shown in Fig. 4. The other formulas on 
DOS and LDOS remain. 

The doping dependence of LDOS has been plotted 
in Fig. 5, where the black solid curve denotes the clean 
case which is equivalent to DOS, and red dashed curve 
denotes the nearest neighbor of the impurity site. It 
has been shown that for all calculated doping, (i) the 
DOS at different doping shows two Van Hove singu- 
larities at corresponding antinodal gap (A„) and max- 
imum gap at the hot spot (A^), reflecting the presence 
of a non-monotonic d^2_y2-waNe gap, (ii) at low ener- 
gies the DOS spectrum also remains (ij,2_j,2-wave like 
and two impurity-induced resonance states in LDOS lo- 
cate symmetrically at positive and negative energies and 
shift towards the low energy with the decreasing dop- 
ing. These features are similar to that shown in the co- 
existing AF SDW and SC order scenario for the some 
small values M (Fig. 3(b) and Fig. 3(c)). However, we no- 
tice that the U-shaped DOS which has been clearly seen 
in the coexisting SDW and SC phase scenario at large 
SDW M can not be obtained in the scenario based on 
dx^^y2-w&Ye with a higher harmonic even at underdoped 
regime (Fig. 5(c)). Our results based on the (i^2_y2-wave 
with a higher harmonic scenario are quite agreement 
with the recent doping dependence of tunneling spec- 



5 



(a) 






A Nr- 








-0.008 



03 (eV) 



FIG. 5: (color online) The LDOS for various doping in the 
clean case and the nearest neighbor of a non-magnetic impu- 
rity site, (a)a; = 0.14, (b)a; = 0.11, and (c)a; = 0.07. A„ and 
Am denote antinodal gap and maximum gap at the hot spot 
respectively. 



trum measurement-, where a non-monotonic gap order 
parameter for the whole doping range studied has been 
shown. Therefore, the non-monotonic gap arising from 
the dx2_y2-weLve with a higher harmonic due to the pe- 
culiar form of the pairing interaction, seems to be more 
reasonable than the physical origin based on coexisting 
SDW and SC phase. 



lyze the DOS and a nonmagnetic impurity effect in the 
electron-doped cuprates superconductor assuming two 
different scenarios: the dx^-y^-'wave superconductivity 
coexisting with AF SDW order versus dx2_y2--wa,ve su- 
perconductivity with a higher harmonic originating from 
the peculiar form of the pairing interaction. We find 
that in both cases there exist two impurity-induced res- 
onance states at low energies. We also find that for the 
intermediate value of the SDW gap, the DOS looks sim- 
ilar to that of the dx2-y2-wave gap with a higher har- 
monic. However, if the SDW gap is sufhcicntly large 
the DOS looks more conventional s-wave like instead of 
the non-monotonic dx2-y2-wave gap behavior. This ob- 
vious difference from that of the c?a,2_j,2-wave gap with a 
higher harmonic even in the underdoped regime, will help 
to differentiate the two above scenarios and make sure 
the physical origin of the unusual non-monotonic gap in 
the electron-doped cuprates. The recent doping depen- 
dence of tunneling spectrum measurement^ provides the 
strong evidence of the non-monotonic gap arising from 
the d2.2_j,2-wave with a higher harmonic, and is qualita- 
tively agreement with our theoretical results. However, 
the earlier point contact tunneling spectrum and pen- 
etration depth measurement show the U-shaped DOS, 
suggesting the coexisting SDW and SC order. Therefore, 
in order to reach the consensus on this topic, our results 
strongly suggest the further doping dependence of tunnel- 
ing spectrum measurement, especially in the underdoped 
regime, is needed to be carried out as soon as possible. 
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